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1.0  INTRODUCTION 


There  is  an  extensive  literature,  spanning  more  than  four 
decades,  describing  investigations  of  spatially  localized 
streams  or  plasma-magnetic  field  entities  (plasmoids)  in 
interaction  with  background  magnetic  field  environments. 
Those  investigations  include  theoretical  analyses,  computer 
simulations,  and  experimental  research  efforts.  The  large 
fraction  of  these  investigations,  however,  have  concerned 
plasma  interactions  characterized  by  MA~1,  where  MA  is 
the  Alfven  mach  number  defined  such  that: 

2  2 

=  (plasma  flow  velocity /Alfven  speed) 

=  (plasma  kinetic  energy  density) /(magnetic  field 
energy  density) 

It  is  well  established  that  in  such  regimes,  the  plasma 
is  either  strongly  decelerated  (and/or  trapped)  by  the 
transverse  magnetic  field  or  propagates  unimpeded  (only 
after  undergoing  charge  polarization  effects)  at  the  reduced 
ExB  drift  velocity.  Very  little  work  has  been  done  in  the 
"super-Alfvenic"  regimes:  MA>>1.  This  is  the  regime  of 
critical  importance  in  the  exoatmospheric  propagation  of 
very  energetic  neutralized  charged  particle  beams.  Recent 
computer  simulations  in  that  regime  (e.g.,  those  of  the 
theory  group  effort  at  Los  Alamos  Scientific  Laboratory 
(LASL))  have  shed  some  light  on  propagation  dynamics  but 
have  been  handicapped  by  the  following  limitations: 


a.  The  lack  of  any  clearly  justifiable  choice  cri¬ 
teria  for  selecting  among  the  many  possibilities 
for  collective  mode  instabilities  in  that  growth 
rates  depend  critically  upon  assumed  initial 
plasma  flow  structures. 

b.  The  computational  complexity  and  time  intensivity 
of  full  three-dimensional  simulations  necessary 
to  treat  correctly  the  boundary  effects  which 
are  inescapable  in  actual  beam  propagation. 

c.  The  limited  ability  to  follow  the  plasrooid  for 
long  (simulation)  propagation  distances  (related 
to  computer  memory  size/computation  time  limita¬ 
tions)  and/or  the  inability  to  choose  rationally 
among  possible  computation  box  boundary  conditions 

d.  An  absence  of  knowledge  of  the  significance  of 
nonlinear  processes  in  either  saturating  linearly 
growing  instabilities  or  aggravating  their  effect. 

Associated  theoretical  efforts  have  almost  always  been 
restricted  to  linearized  treatments  and/or  have  assumed 
(for  tractability)  highly  simplified  geometries. 


It  that  light,  a  carefully  organized  set  of  scaled  simu¬ 
lation  experiments  can  thus  provide  important  new  knowledge 
which  can  support  and  complement  the  theoretical  and  computer 
simulation  efforts.  This  is  the  significance  of  the  effort 
supported  by  the  contract  (F49620-83-C-0091)  now  ended. 

The  experimental  facility  which  was  used  is  located  at 
Morehouse  College,  Atlanta,  Georgia  and  was  made  available 
to  AMAF  Industries,  Inc.  on  a  leased  basis. 
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2.0  STATEMENT  OF  WORK 

The  research  objectives  for  the  effort  as  proposed  were 
defined  by  the  following  tasks: 

1.  Establish  effective  channel  of  communications  with 
existing  computer  simulation  groups  (e.g.,  LASL) . 

2.  Review  published  literature  on  experimental  and  theo¬ 
retical  research  on  instabilities  in  the  plasma  flow  regime 
Mft~l  for  applicable  insights. 

3.  Reconfigure  existing  hypervelocity  plasma  flow  generator 
facility  to  operate  with  high  reproducibility  in  the  regime 
Ma>>1  as  well  as  Mft>l. 

This  will  involve  (as  simple  engineering  tasks) : 

a.  Extension  of  the  vacuum  chamber  in  the  flow  direc¬ 
tion  by  addition  of  a  four-foot  pyrex  glass  pipe 
section. 

b.  Redesign  and  reconstruction  of  the  magnetic  field 
producing  coil  system  (currently  a  Helmholtz 
pair)  to  produce  an  approximately  uniform  field 
over  the  extended  flow  channel. 


c.  Addition  of  a  ballast  chamber,  a  precision  leak 
valve,  a  mechanical  vacuum  pump,  and  a  vacuum 
gauge  to  supply  controllable  quantities  of  selected 
gases  to  the  breech  of  the  gun  (operating  in 
the  gas  "puff"  mode). 

4.  Construct  and  calibrate  fast  response  magnetic,  electric 
and  current  probes  for  flow  diagnostics  (designed  by  well- 
known  techniques,  to  have  response  time  <0.1  ysec) . 

5.  Procure  and  bring  to  operational  status  a  fast  optical 
system  (image  convertor  or  Kerr  optics)  for  macroscopic 
structure  observations. 

6.  Begin  program  of  experimental  parmetric  surveys  (with 
full  diagnostic  monitoring) : 

a.  Plasma  beam  velocity. 

b.  Plasma  beam  density. 

c.  Plasma  beam  particle  mass  (variations  in  composi¬ 

tion)  . 

d.  Propagation  chamber  ambient  pressure. 

e.  Background  magnetic  field  intensity. 


f.  Background  magnetic  field  structure. 


7.  Analyze  results  of  experimental  surveys  as  to: 

a.  Dependence  of  propagation  stability  on  parameter 
values. 

b.  Characteristic  of  instability  structure  (growth 
rates,  spatial  scales,  onset  conditions,  saturation 
levels) . 

8.  Summarize  results  in  form  for  inclusion  in  computer 
simulations  (especially  in  computer  simulation  areas  a. 
and  d.  described  previously) . 

9.  Investigate  experimentally  (guided  by  appropriately 
modified  LASL  computer  code  predictions  as  they  become 
available)  possible  means  for  stabilizing  the  propagating 
beam  with  an  aim  of  achieving  the  highest  energy  in  the 
beam  consistent  with  stability. 
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3.0  REVIEW  OF  PREVIOUS  RESEARCH  RELEVANT  TO  NEUTRAL  PLASMA 
BEAM  PROPAGATION  ACROSS  A  MAGNETIC  FIELD 

This  review  of  previous  research  relevant  to  plasmoid  propa¬ 
gation  is  intended  to  serve  a  number  of  purposes.  Firstly, 
it  presents  results  from  previous  research  which  expounds 
some  of  the  problems  encountered  in  analyzing  plasmoid  propa¬ 
gation.  Secondly,  it  aids  in  the  construction  of  a  founda¬ 
tion  of  basic  concepts  which  must  be  utilized  in  discussing 
plasmoid  propagation.  Lastly,  it  serves  to  motivate  new 
approaches  to  analyze  plasmoid  propagation  with  the  conten¬ 
tion  that  many  of  the  previously  obtained  results  cannot 
be  utilized  in  current  analyses  of  plasmoid  propagation. 

To  begin  the  review,  it  should  be  noted  that  previous  research 
relevant  to  plasmoid  propagation  was  concerned  with  a  wide 
variety  of  applications,  ranging  from  the  interaction  of  the 
solar  wind  with  the  geomagnetic  field  to  the  injection  of 
pulsed,  neutralized  ion  beams  into  devices  designed  to  achieve 
thermonuclear  fusion.  This  wide  variety  of  applicability  has 
presented  problems  to  those  with  present  interests  in  plasmoid 
propagation  in  the  sense  that  it  becomes  difficult  to  extract 
pertinent  information  due  to  the  disparate  parameter  regimes 
involved.  It  thus  becomes  necessary  to  carefully  inspect  pre¬ 
vious  analyses  to  determine  which  parameter  regimes  are  rele¬ 
vant  to  present  interests  in  plasmoid  propagation  (i.e.  mili¬ 
tary  applications) . 
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Another  problem  associated  with  employing  results  from 
previous  research  are  the  limitations  imposed  by  the  highly 
idealized  geometric  configurations  that  were  considered. 

For  example,  Ferraro  ^ ,  concerned  with  the  origin  of 
geomagnetic  storms,  was  able  to  analytically  describe  the 
case  of  a  diamagnetic  effect  -  dominated  semi-infinite 
plasma-field  interaction.  The  limitation  in  this  case 
being  the  absence  of  finite  geometry  (boundary)  effects 
on  the  diamagnetic  current  flow  patterns.  Also,  the  ana¬ 
lysis  performed  by  Ferraro  was  essentially  time- independent, 
which  for  present  applications  is  unsuitable. 

It  can  be  stated  however,  that  Ferraros'  analysis  did  provide 
certain  contributions  that  aided  in  the  understanding  of  plas- 
moid  propagation.  Ferraro,  concerned  primarily  with  the  case 
of  a  plasma  flow  with  sufficient  kinetic  energy  density  to 
exclude  the  ambient  magnetic  field,  demonstrated  that  the 
magnetic  field  decayed  within  a  distance  "c/top"  of  the  surface 
of  the  plasma.  The  quantity  "c/coo",  known  as  the  "collision¬ 
less  skin  depth"  where  "c"  is  the  velocity  of  light  and  "top" 
is  the  electron  plasma  frequency,  is  on  the  order  of  kilometers 
for  the  solar  wind  and  on  the  order  of  a  centimeter  for  high 
density  plasmas  of  interest  for  military  applications. 

Ferraro  has  shown  that  magnetic  field  exclusion  can  be 
expected  to  occur  when  the  square  of  the  Alfven  Mach  number 
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2 

(Ma  ) ,  defined  as  the  ratio  of  plasma  kinetic  energy  density 
to  magnetic  field  energy  density,  is  much  greater  than  one; 

(i.e.  Ma  >>1  ).  For  plasmoids  relevant  to  military  appli¬ 
cations,  the  ratio  (Ma2)  is  on  the  order  of  10*®  and  thus 
serves  as  an  important  parametrical  guideline. 

To  justify  the  importance  of  the  parameter  (Ma2) ,  it  becomes 
necessary  to  trace  the  expected  sequence  of  events  in  the 
flight  of  a  plasmoid  entity.  From  the  aforementioned  kinetic 
energy  density  considerations,  one  is  lead  to  expect  that  the 
initial  phase  of  flight  is  dominated  and  thus  characterized 
by  magnetic  field  exclusion  from  the  interior  of  the  plasmoid. 

The  field  exclusion  is  produced  by  diamagnetic  currents  circu¬ 
lating  in  the  thin  surface  layer  of  the  plasmoid.  These  circu¬ 
lating  diamagnetic  currents  will  eventually  decay  due  to  classi¬ 
cal,  discrete  particle  Spitzer^2^  conductivity  and  also  possi¬ 
bility  due  to  enhanced  anomolous  resistivity  induced  by  non-linear 
instabilities  driven  by  the  diamagnetic  currents.  Thus,  the  decay 
of  these  diamagnetic  currents  permits  partial  diffusion  of  the 
magnetic  field  into  the  plasmoid  after  the  initial  phase  of 
propagation.  An  estimate  for  the  diffusion  time  td  can  be 
obtained  from  the  approximation: 

td  ■  4 ti cl  /C  where  o  is  the  plasmoid  conductivity  and  L 
is  a  scale  length.  For  plasmoids  with  particle  densities  and 
kinetic  energy  densities  relevant  to  military  applications,  and 
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assuming  Spitzec  conductivity,  one  obtains  TD  -  1  micro¬ 
second.  When  considering  the  expected  time-of-flights 
for  military  applications,  and  the  possible  enhancement 
of  anomolous  resistivity,  one  can  conclude  that  the  magnetic 
field  will  partially  permeate  the  plasma  throughout  the 
majority  of  flight.  This  conclusion  has  lead  to  the  develop¬ 
ment  of  several  analytical  models,  attempting  to  describe 
and/or  simulate  plasmoid  propagation  in  parameter  regimes 
relevant  to  military  applications,  which  employ  the  concept 
of  a  polarized  plasma  propagating  across  a  magnetic  field 
through  an  E  x  B  drift  of  plasma  particle  guiding  centers. 

The  validity  of  these  models  is  questionable,  however, 
since  many  ignore  the  strong  plasma  heating  which  may  occur 
during  the  initial  diamagnetic  phase  of  propagation.  It 
is  therefore  the  intention  of  the  next  section  of  this 
report  to  expound  the  limitations  inherent  in  such  models; 
examples  of  which  include: 

3 . 1  MODEL  1:  WILLIAM  PETER  AND  NORMAN  ROSTOKER 
The  objective  of  this  model  is  to  provide  a  theoretical  basis 
for  the  injection  of  pulsed,  neutralized  ion  beams  into 
toroidal  fields  for  the  purpose  of  driving  currents ^  or 
of  supplementary  tokamak  heating.  '  This  model  considers 
only  the  case  of  R  <  a.,  where  R  is  the  plasma  beam  radius 
and  a^  is  the  ion  gyro-radius,  such  that  space-charge  pola¬ 
rization  is  induced  on  the  edges  of  a  bounded  plasma  column 
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as  it  passes  through  a  transverse  magnetic  field.  The  polari 
zation  space  charge  layers  create  a  self-consistent  electric 
field  which  allows  the  plasmoid  to  propagate  across  the 
magnetic  field  by  means  of  the  E  x  B  drift.  The  criteria 
and/or  assumptions  which  characterize  the  conditions  for 
this  motion  to  take  place  are  as  follows: 


3.1.1  The  plasma  beam  must  have  sufficient  kinetic  energy  to 
create  the  polarization  electric  field.  In  analytical  terms. 


1/2  n0Mu02  »_e£ 

8n 

where  nQ  and  UQ  are  the  initial  density  and  velocity  of 
the  plasmoid,  M  is  the  ion  mass,  and  Ey  is  the  polarization 
field.  It  is  further  required  in  this  model  that  the  plasma 
drift  across  the  transverse  magnetic  field  lines  with  its 
initial  forward  velocity  UQ,  which  implies  that; 


2  2 
u)pi  »  0± 
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where  is  the  ion  plasma  frequency,  fti  is  the  ion  gyrofre- 
quency  and  e  is  the  static  plasma  dielectric  constant. 


3.1.2  A  low  beta  ($)  plasma  is  assumed  such  that  the  mag¬ 
netic  field  due  to  induced  plasma  currents  is  negligible. 


Beta,  defined  as  the  ratio  of  plasma  particle  pressure  to 
magnetic  field  pressure,  indicates  the  amount  of  influence 
diamagnetic  effects  have  on  the  plasmoid  propagation. 

3.1.3  The  plasmoid  motion  is  assumed  non-adiabatic,  which 
implies  the  more  general  2-fluid  equations  of  motion  must 
be  utilized. 

3.1.4  All  physical  variables,  such  as  density  and  velocity, 
are  assumed  to  be  functions  of  the  longitudinal  coordinate 
and  time  only. 

3.1.5  Plasma  quasi-neutrality  is  assumed  such  that  the 
ion  density  (n^  equals  the  electron  density  ( n )  (i.e.j 
ni ~  ne  “  where  n  is  a  common  density  called  the  plasma 
density. 

These  assumptions  form  the  foundation  upon  which  the  analytics 
of  this  model  are  based.  Without  belaboring  the  details  of 
these  analytics,  a  review  of  the  quantitative  results  produced 
by  this  model  is  necessary.  This  review  will  enable  compari¬ 
son  of  results  produced  from  other  plasmoid  propagation  models. 

As  a  starting  point,  a  schematic  drawing  of  the  polarization 
of  a  bounded  plasma  beam  in  a  transverse  magnetic  field  Bz  is 
presented  in  Figure  1.  Note  the  coordinate  configuration 


in  which  the  plasmoid  propagates  in  the  x  direction.  Note 
also  that  Ay  represents  the  width  of  the  polarization  charge 
layers.  The  analytical  solutions  for  plasmoid  position  (x) 
and  velocity  (Ux)  derived  for  this  configuration,  neglecting 

O 

terms  of  order  (l/e  ) ,  are: 
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where: 


=  initial  plasmoid  position  for  t  =  0, 
=  initial  plasmoid  velocity  for  t  =  0. 
=  hybrid  gyrof requency  = 

=  static  plasma  dielectric  constant. 
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FIGURE  1.  Schematic  drawing  of  the  polarization  of  a  bounded 
plasma  beam  in  a  transverse  magnetic  field  Bz.  The  resultant 
electric  field  Ex  allows  the  beam  to  propagate  through  the 
magnetic  field  by  means  of  an  EXB  plasma  drift. 


The  oscillatory  behavior  of  the  plasmoid  velocity  in  equation 
[l]  is  due  to  competing  effects  of  the  polarization  electric 
field  and  the  transverse  magnetic  field.  The  frequency  of 
oscillation, \J e  ,  is  the  electron  plasma  frequency  w 

for  e  >>1.  Note  that  the  solutions  are  valid  for  small  values 
of  e  (i.e.;e  t  (M/m)-1-/^)  although  the  quasi-neutrality 
assumption  is  no  longer  valid  in  this  parameter  regime.  For 
a  value  of  e  equal  to  unity,  the  resulting  charge  distribu¬ 
tion  is  unable  to  generate  a  polarization  electric  field  suffi¬ 
cient  to  allow  propagation.  In  this  case,  the  plasmoid  propa¬ 
gates  across  the  transverse  magnetic  field  to  a  scale  length 
equal  to  the  hybrid  gyroradius  aH  and  is  reflected  back  to  the 
injection  plane.  An  expression  for  the  mean  plasma  velocity 
Ux  can  be  obtained  by  averaging  out  the  fast  time  oscillations 


in  equation  £lj .  The  result  is: 
[3]  "x  -  Do  <*  -  ?> 


For  large  values  of  e  (i.e.;t  >>  1),  it  can  be  seen  from  equa¬ 
tion  £lj  that  the  plasmoid  will  propagate  across  the  magnetic 
field  with  a  velocity  nearly  equal  to  its  injection  velocity 
UQ.  In  the  limit  e  >>1,  the  equations  [l]  and  [2]  reduce  to, 
respectively; 


M  ux  =  «o|1  -7  I1  -  cosvC%^-]  1 
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which  is  exactly  the  steady  state  solution  derived  by  Sinelnikov 
and  Rutkevich  in  the  "slow-deceleration  approximation”  to  be 
discussed  in  the  next  section. 

Before  discussing  the  "slow-deceleration  approximation,"  how¬ 
ever,  it  is  worthwhile  to  investigate  quantitatively,  limiting 
values  of  e  and  its  implications  on  the  plasmoid  propagation. 

For  e=l,  we  have  found  that  the  plasmoid  transverses  the  mag¬ 
netic  field  lines  to  scale  length  equal  to  the  hybrid  gyrora- 
dius  (a„)  and  is  then  reflected  back  to  the  injection  plane, 

n 

for  e=10,  the  plasmoid  undergoes  large  oscillations,  and  for 
£  =100,  the  plasmoid  undergoes  smaller  oscillations.  This  is 
depicted  graphically  in  Figure  2.  It  should  be  noted  that  the 
solution  for  the  plasmoid  velocity  Ux  expressed  in  equation 
£4}  is  valid  for  £  >>1,  however,  there  lies  a  further  restriction 
on  the  value  of  e  limiting  the  validity  of  this  model.  Due  to  the 
quasi-neutrality  assumption;  it  has  been  derived  that,  the  domain 
of  validity  of  this  model  is  restricted  to  £  >>(M/m)^^. 

If  this  condition  is  not  satisfied,  it  is  expected  that  longi¬ 
tudinal,  as  opposed  to  transverse,  space  charge  polarization 
effects  will  be  dominant. 
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Plots  of  longitudinal  velocity  Ux/U0  versus  distance  x/a^ 
into  the  field  for  three  different  values  of  the  plasma 
dielectric  constant  e. 

FIGURE  2.  PLASMOID  OSCILLATIONS 

It  is  obvious  that  the  plasma  oscillations  depicted  in 
Figure  2  cannot  be  sustained  indefinitely,  thus,  one  must 
take  into  account  the  expected  damping  of  these  oscillations 
by  introducing  a  collisional  term  into  the  species  (e,i) 
equations  of  motion.  Using  the  Langevin ^  description  for 
the  collision  frequency,  a  term  proportional  to  v (Vj-Ve) 
has  been  introduced,  and  the  result  is  that  the  plasma 
oscillations  damp  out  after  a  few  hybrid  gyroperiods  (ftH  «= 

[ai«e]  V2  >  • 


J 


A 
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Finally,  to  complete  the  review,  a  more  critical  analysis 
of  the  limitations  of  this  model  and  their  implications  is 
necessary.  Firstly,  an  important  limitation  of  this  model 
is  its  inability  to  describe  the  formation  of  the  polariza¬ 
tion  space-charge  layers.  Secondly,  a  more  complete  model 
would  consider  a  truly  bounded  plasma  column,  with  the 
physical  variables,  such  as  density  and  velocity,  being 
functions  of  the  planar  coordinates  perpendicular  to  the 
magnetic  field  direction.  In  actuality,  these  physical 
variables  are  three-dimensional.  Thirdly,  the  extent  to 
which  this  model  considers  longitudinal  space-charge  polari¬ 
zation  effects  is  limited  in  the  sense  that  the  only  results 
yielded  by  this  model  that  are  relevant  to  these  considera¬ 
tions  is  that  the  condition  e  >>  (M/m)^2  corresponds  to 
the  characteristic  distance  for  transverse  charge  separation 
(ai/e )  being  less  than  the  characteristic  distance  for 
longitudinal  charge  separation  (a^ae)l/2.  Lastly,  because 
this  model  considers  low-beta  plasmas  only,  there  also 
exists  an  upper  limit  on  the  values  that  the  static  plasma 
dielectric  constant  e  can  have.  From  the  relation  0  = 

(Uo/C2)e  <<1,  we  have  the  condition  that  e  <<  (Uo/C2)-^. 

Thus  the  limiting  values  of  e  are;  (M/m)1/,2<<  e  <<  (Uo/C2)-1. 

An  explanation  for  the  higher  threshold  e  »  (M/ro)  ^2  required 
for  cross-field  propagation  in  this  model  as  opposed  to 
the  adiabatic  condition  e  >>1  derived  by  Schmidt^,  is 
that  in  this  model,  the  plasma  column  enters  the  applied 
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magnetic  field  from  an  initial  field-free  region;  whereas  the 

condition  e  >> 1  applies  to  plasmas  moving  initially  within  a 
magnetic  field  region. 

3.2  MODEL  2;  K.D,  SINEL'NIKOV  AND  B.N.  PUTKEVICH 
The  objective  of  this  model  is  to  provide  a  theoretical  basis 
for  describing  the  propagation  of  a  bounded  plasma  flux  in 
a  magnetic  field  with  special  emphasis  on  the  formation  of 
transverse  and  longitudinal  polarization  space-charge  layers. 
The  transverse  polarization  space-charge  layers  create  a 
self-consistent  electric  field  (E^)  which  allows  the  plas- 
moid  to  propagate  across  the  magnetic  field  by  means  of  the 
EtxB  drift.  The  longitudinal  polarization  space-charge 
layers  create  a  self-consistent  electric  field  (EL)  which 
may  lead  to  plasmoid  reflection  back  to  the  injection  plane 
depending  on  the  value  of  the  static  plasma  dielectric  con¬ 
stant  e.  The  criteria  and/or  assumptions  which  characterize 
the  conditions  for  this  motion  to  take  place  are  as  follows; 

3.2.1  The  plasma  flux  is  assumed  to  be  bounded  perpendicular 
to  the  magnetic  field  direction  (B)  and  the  plasmoid  velocity 
direction  (V) . 

3. 2. ,2  Plasma  quasi-neutrality  is  assumed  such  that  the  ion 
density  (n^)  equals  the  electron  density  (ne)  (i.e.j  n^-n^n) 
where  n  is  a  common  density  called  the  plasma  density. 
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3.2.3  The  velocities  of  the  electrons  (Ue)  and  ions  (U ^ ) 
are  assumed  to  be  independent  of  the  longitudinal  coordinate 
which  is  parallel  to  the  direction  of  plasmoid  propagation. 

These  assumptions  form  the  foundation  upon  which  the  analytics 
of  this  model  are  based.  Again,  without  belaboring  the  details 
of  these  analytics,  a  review  of  the  quantitative  results  pro¬ 
duced  by  this  model  is  necessary.  This  review  will  enable 
comparison  of  results  produced  from  other  plasmoid  propagation 
models. 

As  a  starting  point,  suppose  that  the  magnetic  field  B  on 

z 

the  x=0  plane  rises  rapidly  from  zero  to  a  constant  value 

Bz  and  remains  constant  for  x>0.  Plasma  bounded  in  the 

y  direction  and  consisting  of  particles  of  mass  m  and  charges 

q  =  +  e  moves  along  the  x  axis  with  velocity  UQ,  and  enters 

the  magnetic  field  at  x=0.  Positive  and  negative  particles 

are  deflected  in  opposite  directions  along  the  y  axis  by  the 

magnetic  field  B  resulting  in  a  transverse  polarization 

z 

> 

electric  field  E^.  Motion  of  the  plasmoid  across  the  mag¬ 
netic  field  results  from  the  Emx  B„  drift. 

In  addition  to  transverse  polarization,  there  exists  in  this 
model,  longitudinal  polarization.  Longitudinal  polarization 
results  from  electrons  being  deflected  by  the  magnetic  field 
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more  than  ions.  This  implies  that  ions  "lead"  electrons  into 
the  transition  layer,  due  to  inertial  effects,  (i.e.;  << 

) ,  resulting  in  a  longitudinal  polarization  space-charge 

> 

layer.  The  resulting  electric  field  E^,  opposite  in  direc¬ 
tion  to  the  plasmoid  velocity,  tends  to  equalize  the  electron 
and  ion  velocities  by  accelerating  the  electrons  and  decelera¬ 
ting  the  ions  such  that  their  penetration  depths  into  the 
magnetic  field  are  equalized,  and  thus,  removing  the  space- 
charge  layer  discontinuity.  Since  plasma-quasineutrality 
is  assumed,  the  relative  longitudinal  displacements  of  elec¬ 
trons  and  ions  is  considered  negligible,  however,  the  resulting 
> 

electric  field  ET  may  be  quite  large  in  magnitude  since 
Me  <<:  Mi*  Also»  the  plasma-quasineutrality  assumption  in 
conjunction  with  the  continuity  equation  implies  that  the 
electron  and  ion  longitudinal  velocities  are  equal.  This 
implication  of  the  quasi-neutrality  assumption  is  justified 
only  if: 


a)  The  electrons  do  not  succeed  in  describing  an 

appreciable  part  of  the  Larmor  orbit  in  a  time 
2 

t  ~  tt/u)0;  where  equals  the  plasma  frequency. 

b)  and,  the  plasma  is  sufficiently  dense. 

If  the  electrons  do  close  their  Larmor  orbits 
in  a  time  t  ~  Vwc»  the  assumption  specified 
in  Section  3.2.3  is  no  longer  valid. 


23 


Now  that  certain  aspects  on  the  formation  of  transverse  and 
longitudinal  polarization  space-charge  layers  have  been 
reviewed,  a  more  detailed  analysis  of  the  "slow-deceleration 
approximation"  can  be  performed.  As  previously  mentioned, 
depending  on  the  value  of  the  plasma  dielectric  constant  e  , 
the  plasmoid  can  be  decelerated  and  then  reflected  back  to 
the  injection  plane  (i.e.;£  ~  1).  For  e  >>  1,  the  analytic 
solution  for  the  plasmoid  velocity  Ux  derived  for  the  "slow- 
deceleration  approximation"  is: 


which  is  identical  to  the  e  >>  1  solution,  derived  by  Peter 
and  Rostoker,  presented  in  Section  3.1.  An  expression  for 
the  mean  plasmoid  velocity  Ux  can  be  obtained  by  averaging 
out  the  fast  oscillations  in  equation  [e] .  The  result  is: 

W  *  uo  [’  -  4-] 

which  implies  that  the  mean  plasmoid  propagation  velocity 
Ux  across  the  magnetic  field  is  not  much  different  that 
the  injection  velocity  UQ  for  large  values  of  e. 

A  graphical  depiction  of  equation  ,  which  is  identical 
to  equation  [4],  is  presented  in  Figure  2;  Section  3.1. 

The  plasma  oscillations  depicted  in  this  figure  cannot  be 
sustained  indefinitely,  thus  one  must  take  into  account 
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the  expected  damping  of  these  oscillations  by  introducing 
a  collisional  term  into  the  species  (e,  i)  equations  of 
motion.  Again,  using  the  Langevin  description  for  the 
collision  frequency,  a  term  proportional  to  vv,  where  v 
is  the  effective  collision  frequency  and  v  is  the  relative 
velocity  of  electrons  and  ions,  has  been  introduced.  The 
result  is: 

a)  if  v<2(0o  damped  oscillations  are  obtained;  and 

b)  if  v>2a>0  an  aperiodic  process  in  which  a  drift 
with  constant  velocity  TJ  =  UD  (1-1/e)  is 
established  in  the  course  of  time. 

To  complete  the  review  of  the  "slow-deceleration  approximation," 
it  should  be  mentioned  that  another  important  objective  of 
the  model  was  to  account  for  the  escape  of  plasma  particles 
from  polarization  layers.  The  escape  of  plasma  particles 
from  polarization  layers,  which  may  be  due  to  the  flow  of 
plasma  along  the  magnetic  lines  of  force,  plays  a  dominant 
role  in  the  issues  which  govern  whether  or  not  practical 
employment  of  plasma  beams  can  be  utilized  for  exoatmospheric 
military  applications.  The  escape  of  particles  from  polariza¬ 
tion  layers  has  an  adverse  effect  on  plasma  beam  directionality 
and  energy  flux  properties  which  are  two  critical  issues  gov¬ 
erning  practical  employment.  Obviously,  the  mechanism  leading 


to  the  escape  of  particles  from  polarization  layers  must  be 
fully  understood. 

To  describe  the  mechanism,  Sinel'nikov  and  Rutkevich  have  as¬ 
sumed  the  plasma  beam  to  be  in  contact  with  electrodes  joined 
through  a  resistance  R.  By  utilizing  the  energy  and  momentum 
balance  equations,  expressions  for  the  polarization  current 
(Ip),  diamagnetic  current  (I\i) ,  and  resistive  current  (IR) 
were  obtained.  From  these  expressions,  they  were  able  to  con¬ 
clude  that  the  plasma  beam  was  decelerated  at  the  entrance  to 
the  space  between  the  electrodes,  which  lead  to  an  inertial 
drift  current  which  compensated  the  current  through  the  resist¬ 
ance  R.  Consequently,  they  were  able  to  determine  that  it  is 
possible  for  the  drift  velocity  Ux  to  be  very  different  from 
the  injection  velocity  UQ,  and  the  energy  of  transverse  mo¬ 
tion  may  be  an  appreciable  fraction  of  the  initial  particle 
energy.  The  appreciable  magnitude  of  transverse  energy  could 
lead  to  plasma  beam  dispersal  which  degrades  the  directionality 
and  energy  flux  properties  of  the  plasma  beam.  If  the  escape 
of  particles  from  polarization  layers  is  due  to  plasma  flow 
along  the  magnetic  lines  of  force,  the  compensating  current 
IR  is  distributed  over  the  entire  thickness  of  the  plasma, 
as  opposed  to  the  narrow  layer  at  the  entrance  to  the  space 
between  electrodes,  and  the  conversion  of  initial  particle 
energy  into  orbital  motion  occurs  along  the  entire  path  of 
motion  along  the  longitudinal  axis.  This  conversion  of  energy 
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into  orbital  motion  again  leads  to  plasma  beam  dispersal  which 
degrades  the  directionality  and  energy  flux  properties  of  the 
beam. 


In  conclusion,  it  is  our  contention  that  the  "slow-decelera- 
tion  approximation"  model  formulated  by  K.D.  Sinel’nikov 
and  B.N.  Rutkevich  is  more  comprehensive  than  the  W.  Peter 
and  N.  Rostoker  model  reviewed  in  Section  3.1,  in  the  sense 
that  it  describes  and/or  simulates  more  physical  effects 
characteristic  of  plasmoid  propagation.  These  physical  effects 
include;  actual  mechanisms  leading  to  plasmoid  propagation 
across  the  magnetic  field,  diamagnetic  response,  the  formation 
of  transverse  and  longitudinal  polarization  space-charge 
layers,  the  escape  of  particles  from  the  polarization  space- 
charge  layers,  and  the  dissipation  of  internally  circulating 


diamagnetic  currents. 
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4.0  W.  R.  SHANAHAN  MODEL 


This  model  combines  a  numerical  and  analytical  study  of 
neutralized  plasma  beam  propagation  across  a  magnetic  field. 
The  model  addresses  physical  questions  of  interest  for  exo- 
atmospheric  military  applications  of  intense  neutralized 
plasma  beams  in  which  attention  is  focussed  on  parameter 
regimes  for  which  the  ratio  of  plasma  kinetic  energy  density 
to  unperturbed  magnetic  field  energy  density  (Alfven  Mach 
Number;  M  )  is  large.  The  review  of  this  model  is  parti- 
tioned  into  three  (3)  sections  as  follows;  1)  Numerical, 

2)  Analytical  and  3)  Final  Results. 

4.1  NUMERICAL  STUDY 

This  section  serves  to  review  Shanahan's  recommended  plan 
of  numerical  study  of  plasmoid  propagation  with  special 
emphasis  on  the  variety  of  computational  tools  employed 
in  the  study.  The  general  scheme  of  the  numerical  study 
envisioned  the  application  of  heuristic,  macroscopic  con¬ 
siderations  and  computations  to  provide  initializing  data 
to  detailed  particle-in-cell  (PIC)  codes  which  were  utilized 
to  simulate  microscopic,  phase-space  dependent  collective 
phenomena.  The  results  of  the  PIC  simulations  were  then 
employed  to  yield  a  refined  description  of  plasmoid  parameters 
which  enter  into  the  macroscopic  treatment  with  the  intention 
of  providing  an  improved  overall  description  of  the  plasmoid 
propagation. 
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To  begin  the  review,  it  should  be  mentioned  that  there 
exists  a  multitude  of  microscopic  phase-space  dependent 
collective  effects  which  can  exert  a  significant  effect 
on  the  macroscopic  behavior  of  the  plasmoid.  These  micro¬ 
scopic  collective  effects  must  be  examined  on  a  time  scale 
roughly  determined  by  the  reciprocal  of  the  electron  plasma 
frequency  wpe,  which  for  typical  parameters  is  approximately 
a  tenth  of  a  nanosecond,  while  total  times  of  flight  are 
expected  to  be  several  milliseconds.  This  disparate  range 
of  relevant  parameters  precluded  the  employment  of  a  single 
analytic  or  numerical  technique  to  perform  the  examination, 
thus,  the  present  plan  of  numerical  study  incorporated 
a  systematic  approach  employing  several  techniques.  The 
conceptual  trend  of  the  systematic  approach,  diagrammed 
in  Figure  3,  was  as  follows: 

4.1.1  The  microscopic  features,  such  as  gross  charge  and 
current  distributions,  of  the  plasmoid  were  determined 
through  calculations  effected  with  a  Magnetohydrodynamic 
(MHD)  fluid  numerical  code.  Classical  transport  coefficients 
were  employed  in  this  initial  phase  of  the  numerical  study. 

4.1.2  Utilizing  the  information  gained  concerning  gross 
charge  and  current  distributions,  initialization  of  PIC 
simulations  was  performed  to  study  microscopic,  phase-space 
dependent  collective  phenomena  which  can  alter  the  transport 
coefficients  appearing  in  the  MHD  fluid  description. 


FIGURE  3.  THE  VARIETY  OF  COMPUTATIONAL  TOOLS 
EMPLOYED  IN  THE  SYSTEMATIC  APPROACH 

4.1.3  The  effect  of  discrete  interactions  of  plasmoid 
particles  with  a  background  neutral  gas  was  then  investigated 
through  employment  of  a  Monte  Carlo  transport  code.  The 
particular  Monte  Carlo  code  employed  in  the  present  numeri¬ 
cal  study  included  the  effect  of  an  externally  prescribed, 
but  not  self-consistent,  electromagnetic  field  on  the  trajec¬ 
tories  of  the  plasmoid  particles.  By  employing  the  Monte 
Carlo  code  with  a  variety  of  initial  conditions,  it  was 
determined  qualitatively,  what  effects  background  interac¬ 
tions  have  on  the  plasmoid  distribution  function.  Data 
acquired  from  the  Monte  Carlo  calculations  was  then  employed 
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to  initialize  the  PIC  simulation  code.  The  objective  of 
performing  PIC  simulations  with  initializing  Monte  Carlo 
data  was  to  evaluate  the  effect  of  the  collective,  restora¬ 
tive  forces  neglected  in  the  Monte  Carlo  calculations. 

4.1.4  Insight  gained  concerning  anomalous,  collectively 
induced  transport  effects  would  then  be  applied  to  yield 
a  refined,  accurate  description  of  the  macroscopic  aspects 
of  neutral  plasma  beam  propagation  across  a  magnetic  field. 

Now  that  the  basic  conceptual  trend  of  the  present  numerical 
study  has  been  reviewed,  it  is  worthwhile  to  describe  the 
various  computational  tools  employed. 

The  principal  numerical  code  employed  in  the  present  study 
of  an  intense  neutralized  plasma  beam  propagating  across  a 
magnetic  field  was  CCUBE,  a  two-dimensional,  relativistic, 
fully  electromagnetic  particle-in-cell  simulation  code. 

CCUBE  can  be  applied  to  study  plasmoid  configurations  in 
any  orthogonal  curvilinear  coordinate  system.  By  employing 
CCUBE  in  the  heuristic  mode,  it  is  meant  that  the  plasmoid 
is  simulated  as  a  whole,  as  opposed  to  specific  region  such 
as  an  ion  sheath  in  which  diamagnetic  currents  may  flow. 

Heuristic,  fluid  calculations  were  effected  with  two  mag¬ 
netohydrodynamic  (MHD)  codes.  A  two-dimensional,  Eulerian 
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code  proved  to  be  the  most  useful  in  the  present  numerical 
study  due  to  the  inherent  features  of  the  code  such  as  the 
inclusion  of  phenomological  resistivity,  energy  loss  due 
to  bremsstrahlung  and  line  radiation,  and  a  moving  reference 
frame  option.  The  second  MHD  code  employed  was  a  two-dimen¬ 
sional  Lagrangian  code  applicable  to  problems  with  fluid- 
vacuum  and  fluid-fluid  interfaces.  The  inclusion  of  phenomo 
logical  resistivity  was  also  a  feature  of  this  code. 

Calculations  concerning  the  interactions  of  plasmoid  parti¬ 
cles  with  background  neutral  gas  atoms  were  effected  with 
the  code  CYLTRAN.  This  code  computes  electron,  ion  and 
photon  transport  consistently  with  the  presence  of  an  exter¬ 
nally  prescribed  electromagnetic  field.  The  code  CLYTRAN 
is  actually  an  embodiment  of  several  Monte  Carlo  techniques. 

Before  completing  this  review  of  Shanahan's  numerical  study 
of  plasmoid  propagation,  it  is  necessary  to  describe  the 
limitations  inherent  in  the  various  computational  tools 
employed.  Firstly,  no  detailed  particle  dynamics  are  in¬ 
cluded  in  the  MHD  fluid  codes,  which  implies  that  in  deter¬ 
mining  the  gross  charge  and  current  distributions,  the 
transport  coefficients  must  be  pre-specif ied.  Arbitrary, 
but  reasonable  values  for  these  coefficients  are  supplied 
by  using  the  collisional  Spitzer  values.  Obviously,  sup¬ 
plying  incorrect  values  would  lead  to  erroneous  results. 
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For  example,  providing  a  value  of  the  conductivity  derived 
from  the  Spitzer  formula  could  lead  to  a  much  smaller  Joule 
heating,  and  accompanying  plasmoid  dispersal,  than  that 
obtained  from  a  more  realistic,  collectively  enhanced  (anoma 
lous)  electrical  resistivity^.  Since  the  principal  quan¬ 
tity  to  be  evaluated  in  Shanahan's  numerical  study  is  the 
collectively  enhanced  resistivity,  which  can  be  determined 
only  by  detailed  microscopic  calculations,  a  cyclic  nature 
is  established  within  the  study  and  is  depicted  in  Figures 
4  and  5.  It  cannot  be  over  emphasized  that  it  is  abso¬ 
lutely  critical  to  obtain  a  refined,  accurate  value  for 
the  collectively  enhanced  resistivity  in  order  to  understand 
the  mechanisms  leading  to  plasmoid  dispersal,  which  is 
one  of  the  major  issues  governing  whether  or  not  practical 
employment  of  intense  neutralized  plasma  beams  for  military 
purposes  is  possible.  Secondly,  restrictions  on  the  appli¬ 
cability  of  MHD  fluid  codes  and  PIC  simulation  codes  arise 
essentially  from  the  two-dimensional  nature  of  the  codes. 
Various  restrictions  include;  the  inability  to  vary  the 
magnetic  field  in  the  third  dimension,  possible  magnetic 
field  exclusion  from  the  path  of  the  plasmoid  is  not  per¬ 
mitted,  current  may  flow  freely  in  the  third  dimension, 
and  consequently  the  magnetic  field  continues  to  increase 
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A  NUMERICAL  STUDY  COMBINING  HEURISTIC, 


MACROSCOPIC  MODELLING  WITH  MICROSCOPIC, 


PARTICLE-IN-CELL  SIMULATIONS 


A  REINFORCING  CYCLIC  FLOW  OF  INFORMATION 


BETWEEN  THE  HEURISTIC  MODELLING  AND  THE 


MICROSCOPIC,  PARTICLE-IN-CELL  SIMULATIONS 


in  the  precursor  region  of  the  plasmoid.  Lastly,  in  the 
PIC  simulation,  a  neutral  plasma  beam  is  injected  into 
a  computational  box  with  specified  boundaries  containing 
a  transverse  magnetic  field.  The  specified  boundaries 
limit  the  plasmoid  time  of  transit  and  in  turn  limit  the 
time  of  simulation.  Also,  the  physical  dimensions  of  the 
plasmoid  must  be  chosen  to  facilitate  fitting  the  plasmoid 
into  the  computational  mesh  without  loss  of  desired  resolu¬ 
tion. 


4.2  ANALYTIC  STUDY 

The  objective  of  the  present  analytic  study  was  to  derive  an 
expression  for  the  dynamic  magnetic  field  response  within 
the  context  of  a  simple  plasmoid  slab  model.  Self-consistent 
coupling  of  this  response  to  the  plasmoid  equations  of  motion 
revealed  the  excitation  of  collective  oscillations.  These 
collective  oscillations,  due  to  initial  charge  separation, 
leads  to  significant  plasmoid  heating  which  causes  plasmoid 
dispersal.  The  dynamic  magnetic  field  response  within  the 
context  of  a  simple  slab  model  was  derived  as  follows: 

4.2.1  Model  plasmoid  as  slab  of  width  L  in  x-direction, 
but  uniform  in  two  transverse  directions. 

4.2.2  Permit  field  quantities  to  vary  in  one  perpendicular 
direction. 
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4.2.3  Derive  equation  for  dynamic  magnetic  field  response. 


Basic  wave  equation 


A2  A  -  -1—  3% 

‘  c2  at2 


where 


Obtain  expression  for  current  (jz)  from  conser¬ 
vation  of  canonical  angular  momentum. 


=*H=r[A°(X°>  -  vx> 


x  *  x  -  v  t 
o  o 


c.  Combine  with  the  basic  wave  equation  to  obtain 
the  desired  equation. 


2  2  2 
_  ko  _  i  a  _  ojp_ 

C*  C2  Y 


]*■■ 


MP  Bq  (X-VQt)  (3) 
C2y 


4.2.4  Solve  equation  (3)  as  an  initial  -value  problem. 
Invert  Laplace  transforms,  and  evaluate  the  asymptotic 
limit  (t  •*  °°)  of  resultant  integrals.  Obtain  expressions 
for  the  dynamic  magnetic  field  response. 


a.  Precursor 
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Field  within  Plasmoid  slab 
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c.  Wake 
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where: 


L  =  Plasmoid  length 
W  =  Plasmoid  width 
VQ=  Initial  plasmoid  velocity 
Kq=  2tt/W 
K2  =  K2  +  wp2/C2 

4.3  FINAL  RESULTS 

The  principal  results  derived  from  the  analytic  study  were 


as  follows 


a.  Wake  Region:  Magnetic  field  reversal  is  exhibited 

immediately  behind  the  plasmoid. 

b.  Slab  Region:  Magnetic  field  exclusion  is  exhibited 

in  the  plasmoid  slab  region. 

c.  Precursor 

Region:  Magnetic  field  compression  is  exhibited 

in  the  front  of  the  plasmoid. 

The  principal  results  derived  from  the  numerical  study 
were  as  follows: 

d.  The  plasmoid  appears  to  propagate  for  short  dis¬ 
tances  (=  2  meters)  without  catastrophic  disrup¬ 
tion. 

e.  A  longitudinal  wave  is  excited  during  this  period 
with  a  wavelength  A  =  10cm.  The  wave  excitation 
within  the  plasmoid  is  due  to  initial  charge 
separation  and  leads  to  significant  heating  in 
the  transverse  direction  parallel  to  the  magnetic 
field. 

f.  Localized  heating  at  the  leading  edge  of  the 
plasma  beam  is  observed  causing  plasmoid  disrup¬ 
tion. 


Before  concluding  this  review  of  the  Shanahan  model  of  plas- 
moid  propagation  across  a  magnetic  field,  it  is  necessary  to 
address  inherent  problems  and/or  anomalies  associated  with 
the  model.  Firstly,  it  is  our  contention  that  mathematical 
errors  have  propagated  throughout  the  analytic  study  in  the 
sense  that  verification  of  the  asymptotic  solutions  for  the 
dynamic  magnetic  field  response  has  not  been  achieved. 
Secondly,  it  is  our  contention  that  the  reference  frame 
transformations  utilized  in  the  analytic  study,  between 
the  laboratory  frame  and  the  plasmoid  rest  frame,  are  not 
proper  in  the  sense  that  the  asymptotic  solutions  for  the 
dynamic  magnetic  field  response  exhibit  anomalous  secularity 
describing  a  linear  growth  in  time.  Therefore,  the  Shanahan 
model,  which  yields  a  wealth  of  information  concerning 
plasmoid  propagation  across  a  magnetic  field  in  various 
parameter  regimes,  must  be  carefully  analyzed  and  corrected 
to  ensure  absolute  validity. 
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5 . 0  STATUS  OF  THE  EXPERIMENTAL  EFFORT 


THE  FACILITY  CONVERSION: 

The  facility  utilized  was  located  in  the  Department  of  Physics 
at  Morehouse  College  in  Atlanta,  Georgia.  It  had  previously 
been  used  in  a  Department  of  Energy  sponsored  research  effort 
in  pulsed,  inductive  direct  energy  conversion  based  on  radial 
expansion  flows.  The  major  components  of  the  system  are: 

A.  The  coaxial  plasma  gun  and  associated  capacitor  bank. 

B.  The  magnetic  field  coil  system  and  associated  capacitor 
bank. 

C.  The  vacuum  chamber  and  associated  pumps. 

D.  Plasma  and  field  diagnostics. 

F.  The  data  acquisition/control  room. 

A  large  area  view  of  the  above-table-top-visible  portion 
of  the  facility  is  provided  in  Figures  6a  and  6b. 

1.  The  Coaxial  Gun  and  Capacitor  Bank  -  The  coaxial 
plasma  gun  is  a  standard  Marshall  gun  configura- 


tion  capable  of  being  either  metallic  foil  or  gas 

"puff"  breech  loaded,  in  the  gas  "puff”  mode,  it 

is  loaded  by  a  thyratron  controlled,  solenoid 

actuated  gas  valve.  The  gun  is  constructed  of 

copper  outer  and  brass  rod  inner  electrodes  with 

an  outer  radius/inner  radius  of  4  and  a  length 

of  10  cm.  The  plasma  composition  options  include 

He,  N,  A,  or  various  metals  (from  foils) .  Details 

of  the  construction  of  the  gun  is  provided  in 

Figures  7a  and  7b.  The  capacitor  bank  for  the 

gun  is  comprised  of  8  capacitors  rated  at  15  uf 

and  20  KV  in  parallel  yielding  a  maximum  storage 

of  24  kilojoules.  The  energy  in  the  bank  in 

switched  by  4  type  GL-7703  ignitrons  in  parallel 

(Figure  8) .  The  capacitor  bank  switched  into  the 

gun  has  a  "ring"  frequency  of  approximately  10  kHz. 

A  typical  storage  oscilloscope  trace  of  the  gun 

voltage  is  shown  in  Figure  9.  Nominal  plasma 

densities  and  flow  velocities  are  respectively 
22  —3 

10  m  and  10  km/sec.  Figure  10  is  an  open 
shutter  photograph  of  the  beam. 

The  Magnetic  Field  Coil  System  and  Capacitor 
Bank  -  The  transverse  magnetic  field  is  produced 


by  a  rectangular  coil  array  extending  most  of 
the  length  of  the  vacuum  chamber.  The  coils 


FIGURE  7a 


FIGURE  7b 


FIGURE  8 


FIGURE  9 


FIGURE  10 


are  visible  in  Figures  6a  and  6b.  The  coils 
consist  of  4  turns  each  and  are  145  cm  long, 

30  cm  wide  and  are  separated  by  25  cm.  The  array 
has  an  effective  inductance  of  approximately 
100  uh.  The  array  is  powered  by  a  capacitor 
bank  comprised  of  4  capacitors  in  parallel  each 
rated  at  60  yf  and  20  KV.  The  bank  is  switched 
by  4  type  GL-7703  ignitrons  in  parallel.  The 
calculated  field  profiles  (normalized  to  the 


field  at  the  array  center)  are  plotted  in  Figures 
11a  and  lib  for  the  orientation:  x,  width;  y, 
length;  z,  separation.  Figure  12  shows  the  array 
location  on  the  vacuum  chamber.  The  capacitor 
bank  "dumps"  into  the  array  with  a  pulse  (criti¬ 
cally  damped)  that  peaks  in  200  psec  and  zeroes 
in  500  ysec.  Figure  13  shows  a  typical  field 
signal  as  picked  up  by  a  magnetic  field  probe 
at  the  array  center.  With  3  KV  on  the  bank, 

the  field  at  the  center  is  measured  to  be  approxi- 

o 

mately  0.05  weber/m  (500  gauss).  The  field 
is  proportional  to  the  bank  charge  voltage  over 
a  large  range  of  initial  charge.  An  estimate 
of  the  fields  required  to  operate  in  the  MA>>1 
regime  is  calculated  in  Appendix  1. 

The  Vacuum  Chamber  and  Pumps  -  The  vacuum  chamber 
is  constructed  out  of  6"  diameter  Corning  pyrex 
sanitary  drain  pipe  sections  with  rubber  gaskets. 
Figure  14  shows  the  overall  dimensions  of  the 
chamber.  The  gun  is  mounted  on  the  axis  of  the 
"T"  end.  The  chamber  is  pumped  by  a  large,  two- 
stage  mechanical  pump  (Welsh  1397)  and  a  small 
oil  diffusion  pump  with  backing  pump  to  the  1- 
10  uHg  range. 
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and  short  time-scale  magnetic  fields  are  measured 
by  magnetic  probes  connected  to  passive  RC  inte¬ 
grators.  Figure  15  shows  such  a  probe  inserted 
through  an  access  port  into  a  crossed-T  section. 

A  typical  calibration  for  such  a  probe  working 
into  a  100  ysec  RC  integrator  is  0.5  mV/gauss 
(5  V/weber/m^) . 

The  plasma  velocity  is  measured  by  "time-of-f light" 
techniques  using  a  bifurcated  light  pipe  and  a 
very  fast  response  photodiode  detector  (Tropel) . 
Figure  16  shows  the  setup  and  Figure  17  shows 
the  location  of  the  velocity  probe  as  well  as 
the  field  probe  on  the  chamber. 

The  density  gradient  structure  of  the  beam  is 
measured  by  a  laser  schliren  system  (Kiefer  & 

Lutz,  1965)  with  a  2  m  travel  path  and  using 
a  5  mwatt  helium-neon  laser  as  a  source.  The 
geometry  of  the  setup  is  shown  in  Figure  18  and 
portions  of  the  actual  system  are  shown  in  Figures 
19a,  19b,  and  19c.  In  this  system  the  voltage 
from  the  photodiode  detector  is  proportional 
to  the  plasma  (local)  density  gradient. 
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FIGURE  19c 


FIGURE  19a 


Figure  19b 
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A  fast  optical  Kerr  cell  shutter  system  has  been 
designed  and  constructed  but  is  not  yet  operative. 
This  system  will  be  used  to  obtain  very  short 
time  interval  photographs  (100  nsec)  which  will 
catch  the  beam  at  early  times  in  its  propagation 
and  allow  for  analysis  of  the  mode  structure 
of  any  dynamic  instabilities. 

The  Data  Acquisition/Control  Room  -  The  control 
room  is  completely  enclosed  by  a  double  screen 
shielded  room.  Figure  20  shows  the  back  of  that 
enclosure  and  the  BNC  data/signal  ports.  The 
control  room  instrumentation  includes  4  Tektronix 
storage  oscilloscopes  (Figure  21a)  for  data  moni¬ 
toring.  Precisely  controllable  trigger  signals 
are  generated  by  the  combination  of  a  Tektronix 
161/162  signal  generator  and  an  Abtronics  Model 
100  4  channel  time  delay  generator  (Figure  21b) . 
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6.0  STATUS  OF  THE  THEORETICAL  EFFORT 

This  section  presents  results  of  the  theoretical  effort 
concerned  with  the  analytic  solution  for  the  dynamic  mag¬ 
netic  field  response  within  the  plasmoid  wake,  precursor 
and  slab  regions.  Two  cases  are  presented  as  follows: 

Case  1:  Includes  inertial  effects  while  excluding 

polarization  and  collisional  effects. 

Case  2:  Includes  inertial  and  collisional  effects 

while  excluding  polarization  effects. 

Case  1  essentially  describes  and/or  simulates  the  conditions 
inherent  in  the  W.R.  Shanahan  formalism  of  plasmoid  propaga¬ 
tion  across  a  magnetic  field.  In  this  case,  the  solution 
for  the  dynamic  magnetic  field  response  within  the  various 
plasmoid  regions  exhibits  a  secular  behavior.  It  is  our 
contention  this  is  an  erroneous  behavior. 

Case  2  addresses  and  attempts  to  resolve  the  secular ity 
issue  by  introducing  a  collisional  term  proportional  to 
the  collision  frequency  (vei)  between  electrons  and  ions. 

The  collision  term  introduces  a  dissipative  mechanism  which 
leads  to  exponential  damping  of  inertial  plasma  oscillations. 
It  is  our  contention  that  in  this  case,  the  secular  behavior 
of  the  dynamic  magnetic  field  response  will  no  longer  be 
erroneously  exhibited. 
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Before  presenting  the  detailed  analytics  of  Case  1  and 
Case  2,  it  should  be  pointed  out  that  our  formalism  attempts 
to  address  several  critical  issues  which  play  a  dominant 
role  in  analyzing  plasmoid  propagation  across  a  magnetic 
field.  These  issues  are;  multi-dimensionality,  collisions 
between  species  (e,i),  transverse  and  longitudinal  polari¬ 
zation,  plasmoid  deceleration  due  to  interactions  with  the 
polarization  fields,  and  the  many  mechanisms  which  lead 
to  the  redirection  of  longitudinal  motion  to  transverse 
motion,  which  in  turn,  leads  to  catastrophic  disruption 
of  the  plasmoid  beam  over  the  time  scale  of  propagation. 

It  now  becomes  necessary  to  point  out  that  the  W.R.  Shanahan 
formalism,  based  upon  the  conservation  of  canonical  angular 
momentum,  which  in  fact  is  not  conserved,  neglects  electro¬ 
static  polarization  fields  in  the  analytic  treatment  while 
simultaneously  incorporating  those  fields  in  the  numerical 
particle-in-cell  simulations.  Likewise,  the  Sinel'nikov 
and  Rutkevich  formalism  considered  to  be  more  comprehensive 
than  the  Peter  and  Rostoker  formalism,  assumes  a  steady 
state  (3/8t-*  0)  in  the  equations  of  motion  for  electrons 
and  ions.  Also,  the  Sinel’nikov  and  Rutkevich  formalism 
does  not  include  a  collisional  exchange  term  in  the  momentum 
balance  equations  of  motion.  The  results  of  our  analytic 
formalism  is  presented  as  follows: 
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Case  1:  Includes  inertial  effects,  excludes  polarization 


and  collisional  effects. 
Simplified  ohms  law: 


£  +  v  x  ^  =  —1 —  31 

uJpc.  3  "t 
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Maxwell's  Equations: 
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Combining  Equations  (1) ,  (2) ,  and  (3)  we  have: 
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Consider  plasmoid  wake  and  precursor  regions: 


(£  .  £.  -  \  Ai  =  O 
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Taking  the  Laplace  transform  of  equation  (8) 


Where: 

We  have: 
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REMARKS  CONCERNING  THE  ANALYTIC  SOLUTION  OF  EQUATION  (13) 

The  analytic  solution  of  equation  (13)  for  the  plasmoid 
wake  and  precursor  regions  describes  electromagnetic  plane 
waves,  originating  within  the  plasmoid  slab,  which  propagate 
into  the  wake  and  precursor  regions  at  the  speed  of  light. 
The  amplitude  (constant)  of  the  plane  waves  is  dependent 
upon  initial  conditions  within  the  plasmoid  slab  region. 


Although  the  actual  analytic  solution  of  equation  (13) 
is  not  presented  herein  due  to  mathematical  uncertainties, 
it  can  be  stated  that  is  our  contention,  the  analytic  solu- 


tion  for  the  dynamic  magnetic  field  response  within  our 
context  (Wake,  Precursor  Regions)  exhibits  an  erroneous 
secular  behavior  if  solved  by  conventional  techniques. 

Note  that  the  analytic  solution  for  the  dynamic  magnetic 
field  response  as  presented  by  W.R.  Shanahan  (Section  4.0) 
also  exhibits  secular  behavior. 

Specifically,  if  the  general  solution  of  the  second  order 
nonhomogenous  differential  equation  (13)  is  represented 
as  the  sum  of  a  particular  solution  Ap  and  a  complementary 
solution  Ajj ,  it  can  be  shown  that  the  general  solution 
of  equation  (13)  for  the  transformed  magnetic  vector  poten- 

tial  (Azwake'  Azprecursor}  wil1  exhibit  secular  growth. 

The  dynamic  magnetic  field  response,  obtained  by  performing 

the  vector  curl  of  the  magnetic  vector  potential,  also 
exhibits  secular  growth.  Note  that  the  nonhomogenous  terms 
(initial  conditions)  of  equation  (13)  arise  from  performing 
the  Laplace  transform  of  equation  (11)  which  contains  a 
second  order  partial  derivative  in  time  and  if  these  non¬ 
homogenous  terms  are  treated  as  driver  terms,  an  erroneous 
solution  will  be  obtained. 

Consider  plasmoid  slab  region:  o*  x'  i.  L 

4  i-v/L*- 


j_  cj.1  A 

cl  r 
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SUU 


6.  x'J  »-v*Vc7 
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(14) 


Taking  the  Laplace  transform  of  equation  (11)  we  have: 


o*: 


s  c* 


(15) 


(16) 


REMARKS  CONCERNING  THE  ANALYTIC  SOLUTION  OF  EQUATION  (16) 

The  analytic  solution  of  equation  (16)  for  the  plasmoid 
slab  region  describes  electromagnetic  plane  waves  bouncing 
back  and  forth  within  the  slab  region  at  a  phase  velocity 
greater  than  the  speed  of  light.  The  group  velocity,  how¬ 
ever,  cannot  exceed  the  speed  of  light.  The  dispersion 
relation  for  electromagnetic  plane  waves  propagating  within 
the  plasmoid  slab  differs  from  the  dispersion  relation  for 

electromagnetic  plane  waves  propagating  in  the  plasmoid 

2 

wake  and  precursor  regions,  by  two  terms  containing  wp 
reminiscent  of  inertial  plasma  oscillations.  These  terms 
in  the  dispersion  relation,  as  can  be  seen  by  referring  to 
equation  (9),  are  due  to  a  driver  term  characteristic  of 
the  plasmoid  current  density  (yQJz) •  It  should  be  pointed 
out  that  the  two  models  of  (Peter  and  Rostoker)  and  Sinel'nikov 
and  Rutkevich)  exhibit  exactly  these  same  inertial  plasma 
oscillations  (Refer  to  Figure  22) . 
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The  solution  of  equation  (8)  for  the  transformed  magnetic 
vector  potential  is  identical  to  that  discussed  in  Case 
1  for  the  wake  and  precursor  regions.  Specifically,  elec¬ 
tromagnetic  plane  waves  propagating  at  the  speed  of  light 
into  the  wake  and  precursor  regions. 

Consider  the  plasmoid  slab  region:  0  i  £.  L 

JT-v.vci 


("  s  '  4Alw\  ^ 


(9) 


Where: 


Equation  (9)  can  be  rewritten  as: 
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Slab 


1  j_  s^60 


\  \fp  &o 


(10) 


Taking  the  Laplace  Transform  of  equation  (6)  and  solving 
for  J 

z 


We  have: 


0'  =  -i-  /  Sk  6. 
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REMARKS  CONCERNING  THE  ANALYTIC  SOLUTION  OF  EQUATION  (11) 
The  analytic  solution  of  equation  (11)  for  the  transformed 
magnetic  vector  potential  is  similar  to  that  discussed  in 
Case  1  for  the  plasmoid  slab  region  except  that  a  dissi¬ 
pative  term  has  been  introduced.  The  effect  of  the  dissi¬ 
pative  term  is  to  damp  the  inertial  plasma  oscillations. 

The  dynamic  magnetic  field  response,  obtained  by  performing 

the  vector  curl  of  the  magnetic  vector  potential  (A  ) , 

z 

does  not  exhibit  secular  behavior  in  this  case.  Again, 
the  actual  analytic  solutions  are  not  presented  herein 
due  to  mathematical  incertainties. 
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COMPUTATION  OF  MAGNETIC  FIELD  LIMIT 


CRITERION:  V/VA  >  10 

PLASMA  FLOW  VELOCITY 
ALFVEN  SPEED 
B/(yo  MN)3* 

B  =  MAGNETIC  FIELD  STRENGTH 

M  =  PLASMA  PARTICLE  MASS 

N  =  PLASMA  PARTICLE  DENSITY 

/.  B  i  V  (y o  MN)VlO 

In  MRS  Units:  y0  ^  10”6 

M  10”26 
N  ^  1022 

V  ^  104 

,  \  B  <  10~2  (102  GAUSS) 


where  V 


V, 
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